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SUMMARY 


Theoretical expressions have been derived by means of linearized 
supersonic-flow theory for the lateral force due to sideslip Cy^, the 

yawing moment due to sideslip C n and the rolling moment due to side- 
slip C 1 for tail arrangements consisting of a vertical triangular 


surface attached to a symmetrical triangular horizontal surface. The 
results are valid, in general, for a range of Mach number for which the 
leading edges of the tail surfaces are swept behind the Mach cone from 
the apex of the arrangement and the trailing edges of the tail surfaces 
are ahead of the Mach lines from the tips . 


A series of design charts are presented which permit rapid estimates 
to be made of the force and moment derivatives. A discussion is also 
included on the application of the expressions for the pressure distri- 
butions determined herein to other plan-form shapes of the tail surfaces 
and possible wing — vertical-tail arrangements. A solution to a two- 
dimensional "mixed type" boundary-value problem which is needed in the 
present analysis but which may also be of interest in other "conical 
flow" analyses is presented in an appendix. 


INTRODUCTION 


The prediction of the stability of complete airplane and missile 
configurations requires a knowledge of the aerodynamic forces and moments 
acting on all the component surfaces of the airframe and the rates of 
change of these forces and moments with the attitude, velocity, and 
acceleration of the associated surfaces. The rates of change of the 
aerodynamic forces and moments when linearly related to the attitudes, 
velocities, and accelerations are commonly called stability derivatives. 
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Theoretical estimates of stability derivatives for a variety of 
wing plan forms with flat-plate cross sections are now available. Infor- 
mation, however, relating to the stability derivatives contributed by 
various nonplanar tail systems is still meager . Most of the available 
derivatives are for configurations composed of low-aspect-ratio surfaces 
(refs. 1 to 3 ) . In reference k, however, sideslip derivatives have been 
presented for tail arrangements for which all the plan-form edges are 
supersonic. In references 1 and 5 approximate estimates of the damping- 
in-roll derivatives for cruciform arrangements with high-aspect-ratio 
surfaces have also been reported. 

The purpose of the present paper is to provide theoretical estimates 
of the lateral force, the rolling moment, and the yawing moment produced 
by the sideslipping motion of a tail arrangement consisting of a triangu- 
lar vertical surface attached to a symmetrical triangular horizontal sur- 
face. The leading edges of the tail surfaces are subsonic; the trailing 
edges, supersonic. Consideration has also been given to the application 
of the results presented herein to other plan-form shapes of the tail 
surfaces and possible wing — vertical-tail combinations. 

The analysis is performed within the framework of linearized 
supersonic-flow theory. Inasmuch as the linearized perturbated flow 
within the Mach cone from the apex of the tail is conical (the arrange- 
ment is a conical body), the analysis reduces to the solution of a sin- 
gular integral equation associated with a two-dimensional mixed type 
boundary- value problem. The solution is obtained by an application of 
the general methods for evaluating these integral equations that have 
been propounded by Muskhelishvili in reference 6 . 

SYMBOLS 


The orientation of the tail arrangement with respect to the X, Y, 
and Z body axes and the positive directions of the velocities, forces, 
and moments are indicated in figure 1. 


X, Y, Z 


body-axes coordinates 


y l» Z 1 


rectangular coordinates in plane parallel to YZ-plane 


V = T} + i£ 

z = x + iy 

0(X,Y,Z) linearized velocity -potential function 

u v w X-, Y-, and Z-components of perturbation velocity, respec- 

tively (v and w are also defined in the v-plane as 
being parallel to the t)- and ((-axes, respectively) 
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z ' = x ' + iy 


u c> v c , W c 

complex velocities, u c = u + iu*, v c = v + iv*, 
and w c = v + iw* 

U*, V*, w* 

harmonic conjugates of the u-, v-, and w- velocities, 
respectively 

V 

free-stream velocity 

M 

free-stream Mach number 


B = ^M 2 - 1 


p 

free-stream density 


free-stream dynamic pressure, ipV^ 

£p 

pressure difference across surface 


pressure coefficient 

a 

angle of attack, radians 

3 

angle of sideslip, radians 

c r 

common root chord of vertical and horizontal tail 

OJ 

fO 

semispan of horizontal tail 

h 

transformed semispan of horizontal tail in y -plane 

\ 

span of vertical tail 

d 

transformed span of vertical tail in v-plane 

x h 

transformed semispan of horizontal tail in z-plane 

x d 

transformed span of vertical tail in z-plane 

e l’ e 2 

arbitrary real constants 

S H 

area of horizontal tail 

Sy 

area of vertical tail 
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7 angle in plane of horizontal tail between a ray through 

origin and X-axis 

7 q angle between leading edge of horizontal tail and X-axis 


t = tan 7 
t Q = tan 7 0 


% 

b 


e 


angle in plane of vertical tail between a ray through 
origin and X-axis 


€ Q angle between leading edge of vertical tail and X-axis 

r = tan e 


tan e 0 


Ay 

2 




ur 

aspect ratio of horizontal tail, = b tan 7 

' c O 

S H 


V 


V 


aspect ratio of vertical tail, —1— = 2 tan 

% ° 


BA 


H 


k = 


2a 2/ tj2 a 2 n 


B Au / B A 


1 - 


H 


B V 


16 / 16 


bV 

16 , 


V 


2a 2 


BA 


1 


2a_ 2 / -d2a 2\ -d2a_2 


BM 


B A H ' B A w BA 


1 - 


H H 


16 / 16 


cn(u/k) 
dn(u/k) 
sn(u/k) 
E, E’ 


Jacobian elliptic functions of argument u and 
modulus k 


complete elliptic integrals of second kind with moduli k 


and \|l - k^, respectively 
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K, K' 


complet e ellip tic integrals of first kind with moduli k 
and ^ 1 - k 2 , respectively 


G = 


BA- 


•H 


E'(k) 


p— K ' (k) + 

k k 


2k 


BV, 

f. B V^ 

1 B 2 A h 2 ba h 

k 1 

r 16 j 

16 “ 4 


Y 

L* 

N 


lateral force, see figure 1 
rolling moment, see figure 1 
yawing moment, see figure 1 


lateral-force coefficient. 




rolling-moment coefficient, 


L' 


ipV 2 b s 
2T TV 


■"n 


yawing -moment coefficient. 


N 




To = 


'*%) 

V^P / p . 


•■0 


c, . *l) 

^ wP / p — >0 


'dc 


n 




^ \<3p / p — *0 

Subscripts: 

H horizontal tail 


V 


vertical tail 
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ANALYSIS 


General Considerations 


The object of the ensuing analysis is to determine the aerodynamic 
pressures and corresponding forces and moments acting on the surfaces 
of the tail arrangement sketched in figure 1 that are produced by the 
sideslipping motion of the tail. The leading edges of the horizontal 
and vertical surfaces are subsonic (within the Mach cone from the apex 
of the system) and the trailing edges are supersonic and at zero angle 
of sweep. It is stipulated that the tail surfaces are of zero camber 
and vanishingly small thickness. It is apparent that this tail con- 
figuration in sideslip attitude is equivalent (by rotation) to a right 
triangular wing at an angle of attack with a triangular end plate or 
fin at zero geometric angle of attack attached to its streamwise edge . 

Such an arrangement is sketched in figure 2, and for convenience this 
orientation of the tail arrangement is considered in the following 
analysis. With the orientation shown in figure 2, the surface approxi- 
mately in the horizontal plane and at a constant geometric angle of attack 
is tentatively called the "wing" and the surface in the vertical plane 
and at zero geometric angle of attack is tentatively called the "fin. 1 

The analysis is based on linearized three-dimensional supersonic- 
flow theory. Specifically, solutions of the linearized three-dimensional 
potential equation 


are sought that satisfy certain boundary conditions associated with the 
wing-fin arrangement. (These boundary conditions are discussed subse- 
quently.) Instead of equation ( 1) , consider the following group of 
equations : 


b2 0xx ' fin ~ fizz ~ 0 


( 1 ) 



(2a) 



(2b) 



(2c) 
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which define "the component-disturbance-velocity fields associated with 
the velocity potential of equation (l) and are more appropriate for the 
ensuing analysis. Once a proper solution to equation (l) or (2) has 
been obtained, the expressions for the lifting pressure may be readily 
determined from the linearized momentum equation 


or 


ipV 2 

2 : 


2 

V dx 


A0(X,Y,Z) 


(3a) 


-A- = | MX,Y,Z) ( 3 b) 

ipV 2 v 


where A0 is the velocity -potential difference across the surface 
and Au is the corresponding longitudinal velocity difference or 
pressure-velocity difference across the surface. Equations (3) are 
consistent with the linearized theory only if the magnitudes of the 
perturbation velocities are equal across the lifting surface. When 
the magnitudes of the perturbation velocities are not equal across the 
^ ^ f t ing surface, equations (3) should contain differences in the squares 
of* the disturbance velocities v and w. The squared terms lead to 
derivatives which are linear functions of a; therefore, these terms 
vanish because our primary interest is the evaluation of the rate of 
change of the aerodynamic forces and moments as a approaches zero. 

Because of the conical geometry of the wing-fin arrangement, the 
following analysis to determine the required solution for the pressure 
employs the concepts of conical-flow theory . This concept implies that 
all disturbance-velocity quantities such as u, v, and w remain con- 
stant along rays emanating from the origin (apex of arrangement) and 
hence become functions of only two independent variables that specify 
the direction of the ray. 

Busemann (ref. 7) initially showed that the assumption of conical 
flow implies mathematically that the disturbance— velocity field within 
the Mach cone of the system is governed by an elliptic differential 
equation, and by a transformation of coordinates this equation reduces 
to the two-dimensional Laplace equation with respect to either the u, 
v, or w perturbation velocity. The problem of obtaining a solution" 
to equation (l) or (2) therefore reduces to one of obtaining a solution 
to Laplace s differential equation in two dimensions subject to certain 
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boundary conditions. These considerations lead naturally to a mode of 
solution using complex-function methods (refs. 8 to 10) and associated 
integral-equation concepts. The following sections present an analysis 
and solution of the wing-fin problem based on these procedures. 


Prescribed Flow Conditions 

A sketch of the wing-fin arrangement showing the body axes used in 
the analysis is presented in figure 3* Denoted also in this figure are 
the prescribed values of the disturbance velocities u, v, and w and 
their spatial derivatives in the plane of the wing and plane of the fin. 
These prescribed values of the velocities and their derivatives are 
determined from a knowledge of the boundary conditions, the symmetry 
conditions, and the equations of irrotationality . 

The boundary conditions are as follows: 

On the Mach cone surface, 

u = v = w = 0 

on the wing surface, 

wfXjYjO 1 ) = aV 

and on the fin surface, 


v 



z) 


= 0 


From symmetry considerations (see ref. 10) it can be shown that in 
the plane of the wing the antisymmetric u- and v-velocities are zero off 
the wing. In the plane of the fin, however, the tangential velocities 
are not zero off the fin because the arrangement lacks symmetry with 
respect to the XZ-plane; in fact, these velocities must be continuous 
across this region. 

The use in the equations of irrotationality of the given boundary 
values of the velocities, together with values of the velocities deter- 
mined from symmetry conditions, produces the additional prescribed 
values of the velocity derivatives denoted in figure 3 and needed in 
the analysis. It is also stipulated that, as the leading edges are 
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approached, the disturbance velocities become locally infinite as the 
7 1 / 2 P ower » that is, the flow around the subsonic leading edges behaves 
in the same manner as the flow around the leading edges of thin flat 
plates in an incompressible flow (see refs. 9 and 10). This stipula- 
tion on the type of edge singularity can be used in order to obtain 
a unique solution to the integral equation of the boundary-value problem 
that is solved subsequently (see appendix A) . 


Transformation of Supersonic Conical Flow to 
Two-Dimensional Incompressible Flow 

The transformation of the supersonic conical— flow equation in u, 
v, or w to the two-dimensional Laplace equation was initially conceived 
by Busemann and expanded in concept and usefulness by many researchers, 
in particular, Lagerstrom, Germain, and Multhopp. Excellent discussions 
of the entire subject are given in the reports by these investigators 
(see refs. 10, 8, and 11, respectively). Only relations pertinent to the 
present analysis are therefore considered herein and the reader is referred 
to the references for proofs and detailed discussions of the relations to 
be presented. 

Figure 4 is a sketch of an arbitrary crossflow plane in the XYZ-space. 
The fact that the body is conical and wholly contained within the Mach 
cone from the apex of the system demands that the u, v, and w dis- 
turbance velocities are homogeneous of zero order and are uniquely defined 
in any crossflow plane j that is, the u-, v-, and w-velocities have the 
functional form 



Now if the X-, Y-, and Z-coordinates are transformed in the following 
manner: 



( 4 ) 
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then equations (2) defining the u-, v-, and w- velocities are transformed 
into elliptic differential equations in the two variables y^ and z^. 

Application of the nonconformal transformation (see ref. ll) 


v 


= T] + iC = 




( 5 ) 


transforms the u-, v-, and w-velocities from the y^z^-plane to the 

v-plane in which each of the disturbance velocities satisfies Laplace's 
equation; that is, in the v-plane, 

sp u = 0 

X^v = 0 
^w = 0 

o /^2 ^2 \ 

where v is the Laplacian operator [ + — - . The effect of the 

W d£ 2 / 

transformation expressed by equation ( 5 ) is to deform the doubly con- 
nected annular region between the Mach cone and the body in the y^-plane 

(see fig. 4 ) into a doubly connected open-slit region in the v-plane (see 
fig. 5 ). The continuous cut (l,oo,-l) in the v-plane corresponds to the 
Mach circle in the y^^plane. If the v-plane is considered as composed 

of two sheets, then the region external to the Mach circle in the 
y-^-plane transforms into the lower sheet of the v-surface and is con- 
nected to the upper sheet through the branch cut (l,co,_l), that is, the 
transformed Mach circle. The transformation does not distort the hori- 
zontal and vertical axes, and hence the shape of the wing-fin arrange- 
ment is preserved (except for scale) in the v-plane. 

It should be pointed out that, since the wing-fin contour is trans- 
formed to one sheet (the upper sheet) of the double-sheeted v-surface, 
the correspondence between the velocities in the crossflow plane of the 
original XYZ-system within the Mach circle and the upper sheet of the 
v-surface is 1:1. Furthermore, since the transformation is conformal in 
the neighborhood of the q- and £-axes, the prescribed values of the 
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velocities and their spatial derivatives which are constant along the 
wing plane and fin plane in the original space remain constant along 
the t]- and £-axes in the v-plane. 

If the complex velocities u c = u + iu*, v c = v + iv*, and 

w c = w + iw* are considered in the v-plane, then from the analysis of 
Hayes and Multhopp (refs. 9 and 11) the ensuing compatability relations 
(which take the place of the equations of continuity and irrotationality) 
provide the necessary relations in the v-plane for attempted solutions: 


dv 


c 



(6a) 


dw 


c 



(6b) 


In terms of the real parts of u c , v c , and w c , the relations of equa- 
tions (6) are as follows (see ref. ll) : 

For £ = 0 and -1 < tj < 1, 


du _ dv 
dt] df] 


B — - H dw 

" " fiTlp ^ 

B — - T) dw 

dr < = {TT^s 

dv _ 1 dw 


( 7 ) 
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and for q = 0 and £ ^ 0 


dw _ g ) 1 + du 




dv -g 1 du 

d^ £ 8 t] 


(8) 



dv -o 1 du 


By using these relations, general solutions for the u, v, and w dis- 
turbance velocities can be determined, provided a solution exists for 
any one of the velocities. 


The evaluation of the u perturbation velocity in the v-plane 
takes on its simplest form when the complex sidewash velocity v Q along 

the wing-fin contour is initially determined and the u-velocity is then 
derived from the sidewash velocity by using equations (6) or (7). The 
expression for the sidewash velocity along the contour has been evaluated 
in appendix A by obtaining a solution to the integral equation defining 
the sidewash in terms of its prescribed boundary values. 

In appendix A the expression for v c is initially derived in a 
z-plane which is obtained conformally from the v-plane by the following 
transformation: 


where the plus sign is valid for x > 0 and the minus sign for x < 0. 
Figure 6 is a sketch of the z-plane and shows that the wing-fin contour 
of the v-plane becomes a slot along the real axis (the x-axis) of the 
z-plane. The details of the transformation from the v- to the z-plane 
are given in appendix A. In the z-plane the expressions for the u-, 
v— , and w- velocities can be expressed much more compactly and simply 


Evaluation of the u Perturbation Velocity 
(Pressure Velocity) in the v-Plane 
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than in the v-plane and for this reason many of the following solutions 
for the velocities are given in this plane . The expressions for the 
real and imaginary parts of the complex sidewash velocity along the 
transformed contour in the z-plane as derived in appendix A are as 
follows: 


( 9 ) 

(10) 

where e-j^ and e 2 are arbitrary real constants. 

The u-velocity is immediately determined from the following rela- 
tions in the v-plane that can be determined from equations (6): 


v (x ,on . K ♦ - 


X - X. 


X* - X 


(x h ^ X < X*) 


v *(x,0 + ) - (e, ♦ 


(— X h ^ X ^ x h ) 


du dv 


(- 1 < T ] < 1 } 5 = 0 ) ( 11 ) 


B — - £ dv* 

d£ d 5 

In the z-plane these relations become 


(n = 0; 5^0) (12) 


B 


du(x,0+) 

dx 


= - \Jl + h 2 \jx 2 - X h 2 — ( ^° +) (x h g X g x d ) 

ox ' 


(13) 


du(x,0t ) 
dx 


- _\li j. \[72 ~2 dv*(x,0+) , ^ ^ \ 

\ji + n yx h - x ^_x h g x g Oj 


> 


(ii*) 
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Equation (13) relates the u- and v- velocities along the transformed 
wing contour in the z-planej equation (l4) relates the v*~ and 
u _ velocities along the tr an s f orme d fin contour in the z— plane • 
Integration of equations (13) and (l4), with equations (9) and (10) 
for v and v* taken into account, produces the following solution 
for the u- velocity: 


u(x,0 + ) = 


M + h c 


B 




.( x d - \) 


X - X 


+ e. 


h 


x h + x 



x h) 


+ 


log c 


x d - x h - 2x + 2^|(x h + X)(x - x d ) 


x d + x h 



(“ x h ^ x ^ x d) (15) 


The relationship between e 1 and e 2 is easily determined by using 

equation (15) in the following manner: In the v-plane u(0“,0 + ) is 

known to be zero from the physical boundary conditions. Hence in the 
z-plane u(x,0 + ) is zero at x = -x^. In equation (15) for the 

u-velocity , if x is set equal to -x h and u^Cr*") is set equal 
to zero, the following relation between e-j_ and e 2 is obtained: 


e 


1 



(16) 


The same relation between and e 2 is also obtained from the con- 
dition that the circulation along a closed path enclosing the wing-fin 
contour must be zero (see ref. ll) . The substitution of e^ in terms 

of e 2 , into equation (15) yields 
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u(x,0 + ) = 


\ 


— ■ h e 2 (1 + Xh + x 


B 




X / \ x d - X 


^ x ^ x d) (17) 


The constant e 2 has been evaluated in appendix B and may be simply 
expressed as 


eo = 


Va 1 


1 + h c 


where 


G = 


x h K'(k) + 


E'(k) 


2k 


x d " x h 


_ 1 - va - \lh - x h 2 )( i - 


x d - x v 


In the equation for G, K'(k) i s the complete elliptic integral of the 
first kind with modulus \Jl - k 2 , and E’( k) is the complete elliptic 

integral of the second kind with modulus jl - k 2 . The parameter k 
necessary for the determination of G may be expressed in terms of the 
geometric characteristics of the tail system (see appendix B) . These 
relationships are illustrated in figure 7. 

The u-velocity along the contour in the z-plane is now uniquely 
given as J 


u(x,0 + ) = (1 + — id 


BG 


*h 


+ x 


X* - X 


(18) 


• Ttle ex P res ® ions ^ or "the u-velocity along the wing and fin contours 
in the v-plane (see fig. 5) are easily obtained from equation (l8) through 
the use of the transformation 
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z 


= x + iy = ± 


2 , ,2 

v + h 

1 + h 2 


( 19 ) 


where, as previously stated, the plus sign is valid for x > 0 and the 

sign for x < 0. In the v— plane the expressions for the u— velocity 

along the contours become: 

Along the wing (0 ^ T) S d) , 


(l,0*) =±^[x + 


h 


h + \ T| 2 + h 2 


EG y \Jt] 2 + h 2 ^d 2 + h 2 - \] T ) 2 + h 2 

along the upper half of the fin (0 = £ =h), 


( 20 ) 



A + , h 

V 


h + \ h 2 - £ 2 


d 2 + h 2 7 


\|h 2 - 'i‘ 


(21) 


and along the lower half of the fin (-h S ( = 0)> 


u(S,0*) 



h ._A 

h + )|h 2 - ^ 

<|h2 - S 2 /\ 

'Jd 2 + h 2 ± \Jh 2 - £ 2 


( 22 ) 


The equation for u on the lower half of the fin (eq. (22)) has been 
obtained from equation (2l) for the upper half of the fin by applying 
the known condition of antisymmetry of the u-velocity with respect to 
the q-axis . 


The expressions for the pressure distributions for the wing and 
fin follow imm ediately from equations (20) to (22) through the use of 
the following linearized pressure-velocity relation: 


p = -pVu 


( 23 ) 
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where p represents the disturbance pressure. The use of equation (23) 
yields the following expressions for the pressure coefficient for the 
wing and the fin: 

On the wing (0 ^ t) £ d) , 

tip p(t), 0 - ) - p(t),O f ) 


|pv 2 


ipv 2 


= |fu(^0-) + ud,,^ 


h + 


T) 2 + h 2 


^/l + *L_J 

BG \ \^2V h2 / |[ \jd 2 + h 2 - + h 2 


(24) 


on the upper half of the fin (0 ^ ^ ^ h), 


tip • p(£,o“) - p(^0 + ) 


|pv 2 


i pv2 


= f £U,0 + ) - u(C,0-[) 


2a 

/1 + h v 

oj 

1 

+ 

BG 


\ ^d 2 + h 2 - \jh 2 - t 2 


/1 - h v 

h - fh 2 - t , 2 

V k 2 - m 

^d 2 + h 2 + 'jh 2 - t 2 


( 25 ) 
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and on the lower half of the fin (-h ^ C ^ 0) 



h - ^h 2 - C 2 


(Jd 2 + h 2 + \| h 2 - t 2 


2 



h +\|h 2 - £ 2 


(26) 


^d 2 + h 2 - ^h 2 - t 2 


Evaluation of Pressure Relations for the Vertical- and 
Horizontal-Tail Combination in the Original Space 

The results of the preceding analysis (in particular, eqs . (24). 
to (26)) allow the evaluation of the pressure coefficients for the wing- 
fin arrangement in the original XYZ-space. It is convenient, however, 
at this point to return to the orientation of the arrangement that is 
shown in figure 1 so that the wing now becomes the vertical tail and 
the fin the horizontal tail. It should be recalled that for convenience 
in the preceding analysis the tail arrangement of figure 1 was rotated 
to the position shown in figure 2 and the vertical tail (the lifting 
surface) was tentatively called the wing and the horizontal tail. (the 
end plate) was tentatively called the fin. With this point in mind the 
following expressions for the pressure coefficient apply to the vertical, 
and horizontal-tail combination orientated as shown in figure 1 and obey 
the sign conventions indicated therein. The expressions for the pres- 
sure coefficients in the original XYZ-space are obtained from equa- 
tions (24), (25), and (26) by the following transformation of the T},t 

variables (fig. 5) to the |,| variables by using equations (4) and (5) 
In the plane of the vertical tail (or wing), 



(27) 


d = B — = B tan e Q = Br Q 
c r 


BA y 

~ 2 ~ 


(28) 


where r is the slope of a ray from the apex, by is the span, 
and Ay is the aspect ratio of the vertical tail. 

In the plane of the horizontal tail (or fin), 


is the semiapex angle, 


C = 


B - 

X 


B tan 7 


Bt 


1 _ B 2 ^ 2 \Jl - B 2 tan 2 7 \l _ B 2 t 2 


( 29 ) 


h = 



B tan 

' 0 

Bt 0 

Att 

B f 


- -w 

^1 - B^tan^7 Q 

■I 1 - B V 

|l - B 2 

a h 2 

15" 


( 30 ) 


where t is the slope of a ray from the apex, bg/2 is the semispan, is the semivertex 

angle, and A H is the aspect ratio of the horizontal tail. The substitution of equations (27) 

to (30) into the pressure-coefficient equations (2b) to (26) results in the following relations 
for the pressure coefficient for the original tail arrangement in the XYZ-space of figure 1: 

For the vertical tail, 



Bt, 


B 2 t r 


B 2 r 2 + 


B 2 t Q 2 


1 - B^t, 



Bt, 


\|1 B 2 t 0 2 


B 2 r 2 + 


B 2 t Q 2 
1 - B 2 t Q 2 


B 2 r Q 2 + 


B 2 t Q 2 
- B 2 t Q 2 


( 31 ) 


B 2 r 2 + 


B 2 t 0 2 
1 - B 2 t Q 2 


for the left half of the horizontal tail. 


ro 

o 



(32) 


hp 
I p v 2 J 

2 P / H 


I 





(33) 
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The parameter G appearing in these expressions 
figure 8 wherein the variation of — 


16 

B V 


may he determined from 


- 1 G with BAy is 


presented for different values of 
mations of G from figure 8, the 
BAy for various values of BA H 


BAtt. In order to facilitate esti- 


variation of 


B 2 A V 2 


16 


IB 2 A H 2 


- 1 with 


is presented in figure 9. 


An idea of the spanwise variation in pressure and resultant pressure 
for a typical vertical- and horizontal-tail combination may be obtained 
from figure 10. These results were obtained by use of equations (31) 
to (33). 


It is of interest to know the expressions for the pressure coeffi- 
cients for limiting cases of the tail arrangement obtained, for example, 
by setting the horizontal-tail span equal to zero or letting all the 
leading edges coincide with the Mach cone traces in the plane of the 
tail surfaces (sonic leading edges). Such expressions are obtained from 
the general expressions - equations (31) to (33) - for the pressure coef- 
ficients, provided the limiting values of the parameter G are avail- 
able. The values of G for a number of limiting cases of the tail 
arrangement have been evaluated. These values of G, together with the 
corresponding expressions for the pressure coefficients, are presented 
in table I. It should be noted that results presented in table I for 
various limiting arrangements are identical to previously published 
results (refs. 4, 12, 13, and 14) that are available for some of these 
limiting arrangements . 


Evaluation of Forces and Moments 

The resultant aerodynamic forces and moments acting on the tail 
arrangement in a sideslip attitude are easily obtained by use of the 
pressure-coefficient expressions given by equations (31) to (33). The 
expressions for the resultant forces and moments to be determined obey 
the standard sign convention indicated in figure 1. 
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Lateral-force derivative Cy^.- For the 
tive sideslip attitude, the lateral force is 


tail arrangement in a posi- 
defined as 


Y = - 



dx dz 


( 5*0 


where ( — ) is the pressure coefficient for the vertical tail and Sy 

\< 1/V 

represents the vertical-tail area. The minus sign is prefixed to the 
integral so that the lateral force obeys the sign convention indicated 
in figure 1. It is clear that the horizontal tail gives a null contri- 
bution to the lateral force. The lateral-force derivative is defined as 



d / Y \ 

^3 \qSy J 


0— >o 





dx dz 


0 


(35) 


The substitution into equation (35) of "the expression for obtained 

from equation (31) yields, after integration, the following expression 
for C Yp : 




ij-n 


B c 


( 3BA H 

bV bV \ 

\ji6 - b 2 a h 2 1 

^ 16 - B 2 A h 2 J 


(36) 


As stated previously, estimates of the factor G can be obtained from 
figures 8 and 9. The variation of BCy^ with BA y for various values 

of BA H is presented in figure 11. 
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Rolling-moment derivative Cj^.- The rolling moment produced by 
a positive sideslip attitude of the tail is defined as 



where Sg represents the horizontal-tail area. Note that the horizontal 
tail contributes a positive moment to the total rolling moment. The 
rolling -moment derivative is defined as 




a / l' 


^3 \qSybyi 


— p— >0 


C l R =1 


dp 


Syb v JJ s v W/V 


(?: 




z dx dz + 


Syby JJ s H \ 


y dx dy 


0 — >0 


% = (%)v + (%) 


P/H 


(37) 


Substituting into equation (37) the expressions for ( — | and ( —) 

\ 9 - / V \ ^ /H 

obtained from equations ( 31 ), ( 32 ) > and ( 33 ) and then performing the 
necessary integrations yields 


C Gz p) v 


3Gb5a v 2 


7BA h 


+ 3' 


Jl6 - B ^ H 2 


A. 2 *2 a 2 


k 16 - b 2 a h 2 


2BA W 


J16 - b 2 a h 2 


b 2 a v 2 


, 2 a 2 


B*Ai 


BA„ 


* 16 - B V \ji6 - bV/ 


2\ll 4 


i 2 a v 2 b 2 a h 2 


ba r 


16 - b 2 a H 2 ^16 - B 2 A h 2 I 


I + sin - 1 


B 2 A v 2 B 2 A h 2 


3 BA h 


16 - B 2 A h 2 \Jl6 - B2AJJ2 


B 2 Ay 2 B^Ajj 2 BAjj 


r 1- * 


16 - B A h ^1x6 - B 2 A h 2 j 


m 
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and 




( 59 ) 


Variations of B ( c Zp) v and B ( C Zp)n with BA v and are P resen ' fc e ( 3- 

in figures 12 and 13 , respectively. 


MCA TN 3071 


25 


Yawing-moment derivative C n ^ .~ The yawing moment about an axis 
through the apex of the arrangement is given by 


N = 



m 


(the horizontal tail offers no contribution to the yawing moment) \ 
N may also be written as 


N = xY 


where x is the x-coordinate of the point at which Y acts (often 
referred to as the center of pressure) » An examination of the equation 

for (31)) shows that it is conical in form; hence, the 

point x will be 2. c r for a triangular vertical tail. This statement 
may be expressed in equation form as 



The yawing-moment derivative may now be simply expressed as follows: 



- — Cy 

3A V ^3 


or, by utilizing equation (36), 





l6rt 


3BA H 


2a. 2 


BM 


2a 2 


\Jl6 - B 2 Ah 2 1 


16 - B 2 A h 2 


(^3) 
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For convenience , the variation of 

of BAg is presented in figure Ik. 

other position of its moment axis, 
needed: 


with BA y for a range of values 


ULp vv-LOii 

In order to determine 


C n p for some 


the following transfer formula is 


Cn 


P 



(MO 


where l is the distance parallel to the x-axis between the new position 
of the moment axis and the apex of the tail. Considering equation (4l) , 
equation (44) may also be written as 



(45) 


DISCUSSION OF RESULTS 


In the foregoing section general expressions for the sideslip deriva- 
tives C Yp (eq. ( 56 )), C lp (eqs . ( 57 ) to ( 59 )), and C np (eq. (43)) 

were derived. It is of interest to present results for some limiting 
cases of these expressions. Table II contains analytical expressions 
for these derivatives for the same limiting cases of the tail arrange- 
ment that were considered for the pressure coefficients presented in 
table I. For the limiting case BAy = 2, that is, for a vertical 

tail with a sonic leading edge (the Mach line coincides with the leading 
edge), the variation of BCyp, BCjp, and C np with BAg is presented 

in figure 15 . In a similar manner the variation of these derivatives 

with BAy for BAg = 4 (the leading edges of the horizontal tail are sonic) 

is presented in figure l6. From figure 14 the result that the yawing- 
moment-coefficient derivative of the vertical tail C np is finite for 

BAy =0 is at first sight surprising. One would expect all derivatives 

to vanish when the disturbance surface (the vertical tail) vanished; 
however, it is to be recalled that the derivative C np is made non- 

dimensional with respect to Syby and, although the yawing moment N 
vanishes as Sy > 0, the ratio N/Syby which is in essence the coeffi- 
cient remains finite since by >-0 as Sy >• 0 and Syby gives an 

infinitesimal of the same order as N when Sy > 0. 
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An indication of the variation of Cy p , C Zp , and C n with Mach 

number for fixed aspect ratios of the vertical and horizontal tails can 
be obtained from figures 17, 18, and 19, respectively. The results pre. 
sented in these figures for supersonic leading edges (the Mach lines 
behind the leading edges) were obtained from reference k. 

Although the results of the analysis presented herein are for a 
tail arrangement consisting of triangular horizontal and vertical sur- 
faces, the results could be easily extended so as to include tails with 
sweptback or sweptforward trailing edges but of zero taper (as sketched 
in fig. 20), provided the angle of sweep of the trailing edges is less 
than the sweep angle of the Mach lines from the apex of the trailing 
edges (supersonic trailing edges). This limitation follows of course 
from the classic results of linearized flow theory that disturbances 
propagate only downstream and the pressure coefficients presented 
herein for the triangular surfaces would also be valid for sweptback 
or sweptforward trailing-edge surfaces with the stated limitation. 


It is of interest to note that the expressions for the pressures 
presented in table I for the case where the leading edges of the hori- 
zontal tail are sonic are also valid for the wing— vertical-tail 
arrangements sketched in figure 21, provided of course it is realized 
that the wing must be at zero angle of attack and that the Mach lines 
from the apex of the vertical tail that are in the plane of the wing 
must intersect the trailing edge of the wing. 


CONCLUDING REMARKS 


Application of linearized theory has enabled an evaluation of the 
lateral force due to sideslip Cyn , the rolling moment due to side- 


slip 


c v 


and the yawing moment due to sideslip C n for a tail 


arrangement consisting of a vertical triangular surface attached to a 
symmetrical horizontal triangular surface. A series of design charts 
have been prepared which permit rapid estimates to be made of all the 
derivatives for a range of Mach number for which the leading edges are 
subsonic and the trailing edges are supersonic. 


The expressions for the pressure coefficients determined, in addi- 
tion to being valid for the plan forms considered herein, are also valid 
tail arrangements with sweptback or sweptforward trailing edges but 
of zero taper, provided the angle of sweep of the trailing edges is less 
than the sweep angle of the Mach lines from the apex of the trailing 
edges . 
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The expressions for the pressures also apply without change to cer- 
tain wing— -vertical-tail arrangements, provided the wing is at zero angle 
of attack and the Mach lines from the apex of the vertical tail that are 
in the plane of the wing intersect the trailing edge of the wing. 

The solution presented for the two-dimensional "mixed type" boundary- 
value problem and needed in the present analysis is fundamental for prob- 
lems involving vertical- and horizontal-tail combinations and is suitable 
for application to more complex systems than those considered herein. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., October 16, 1953* 
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APPENDIX A 


EVALUATION OF THE LINEARIZED DISTURBANCE VELOCITY 
ALONG THE CONTOUR OF THE WING-FIN 
ARRANGEMENT IN THE v-PLANE 


The evaluation of the tangential disturbance velocity along the 
wing-fin contour in the v-plane can be performed by an application of 
techniques discussed by Muskhelishvili in his book on singular integral 
equations (ref. 6). 

First, the v-plane is conformally transformed into the z -plane 
by use of the transformation 


v 2 + h 2 
1 + h 2 


(Al) 


where the plus sign is valid for x > 0 and the minus sign is valid for 

x < 0. In the z-plane (see fig. 6) the image of the wing-fin contour is 

a slot along the real axis (the x-axis). The branch cut (l,oo,-i) in the 
v-plane representing the original Mach circle retains itself along the 
real axis in the z-plane . 

Since the z-plane is a conformal map of the v-plane, the complex 
velocities u c = u + iu*, v c = v + iv*, and w c = w + iv* which satisfy 

the Laplace equation in the v— plane also satisfy the Laplace equation in 

the z-plane. The given boundary values of u, v, and w and their 
spatial derivatives which are constants along the real and imaginary axes 
in the v-plane therefore remain constant in the z-plane along the real axis. 

will be found expedient to determine the variation of v— velocity 
(the real part of v c ) along the wing -fin image (from x = -x^ to 

x = x d) * Inasmuch as this lateral velocity is antisymmetric with respect 

to the q-axis in the v-plane and the wing-fin contour is symmetric with 
respect to the t]— axis, it is necessary to consider a solution for v in 
only the upper half of the v— plane and therefore in only the upper half 
of the z-plane. Figure 6(b) is a sketch of the upper half of the z-plane 
and. shows the appropriate boundary values of v and v* along the real 
axis . From the boundary condition that w = cxV on the wing, it is 
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possible to determine that v* = 0 (or constant) along the wing segment 
x h ^ x g x d in the following manner. By using equations (6a) and (6b), 

the relation in the v-plane between v c and w Q is found to be 


dv c = i - - 1 - - dw c (A2) 

\|l - v 2 

Along the wing axis v = t), and by using equation (A2) can be 

expressed as 

dv* = — dw ( -1 ^ T) ^ l) (A 3 ) 

Now w = Constant along the wing segment 0 ^ r] ^ d or dw = 0; hence, 
v* = Constant. This constant may be set equal to zero (see ref. 6). 

This condition v* = 0 remains valid along the image of the wing segment 
in the z-plane . 

There is one factor in the transformation of the wing-fin contour 
from the v-plane to the z-plane that requires particular attention. 

Note that the point ih representing the leading edge of the fin in 
the v-plane transforms into the origin in the z-plane and that the two 
faces of the fin spread out along the real axis of the z-plane, one on 
each side of the origin, and form a continuous segment (-Xh>°> x h) . Now 

if the origin is approached from above along the line x = 0, the 
v-velocity will become infinite in the immediate neighborhood of the 
origin even though the origin does not represent a free edge such as 
x^ where singularities in the boundary functions are usually allowed. 

This discontinuity in v follows directly from the behavior of the 
v-velocity in the neighborhood of the fin leading edge ih in the 
v-plane. In this plane as the fin edge ih is approached from above 
along i£, the v-velocity becomes a singularity of order -l/2; hence, 
in the z-plane as the origin (the image of the point ih) is approached 
from above along iy, the v-velocity becomes a singularity of order -1. 
The fact that the singularity in the v-velocity near the origin in the 
z-plane is of -1 order arises from the double-valuedness of the trans- 
formation from the v-plane to the z-plane. For the subsequent analysis 
the point to remember from the present discussion is that the v-velocity 
function does not continuously approach its prescribed values for all 
points of the open segment -x h < x < x h but becomes discontinuous as 

the origin is approached along x = 0. Along the remaining segments 
defining the wing-fin image in the z-plane, the v-velocity can be shown 
to approach its boundary values continuously except at free edges where 
singularities are of course allowed. 
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As indicated previously, the solution for v c is based upon the 

analysis of reference 6. In order to associate the present analysis 
directly with the general treatise on boundary-value problems in refer- 
ence 6, one additional ■transformation is needed^ this transformation is 


2z - (x h + xj 


(A4) 


In f igure 22 the transformed wing-fin contour in the z ' -plane is sketched 
and the corresponding boundary values of v along L' and v* along 
L' are noted. It should be noted that only the upper half of the 
z -plane is being considered with boundary conditions prescribed along 
the axis of reals . The problem at hand of determining v c may be named 

the two-dimensional mixed- type boundary-value problem for the half plane. 

It is well-known (see refs. 6, 8, and ll) that the solution to this 
type of boundary-value problem essentially involves the solution of the 
singular integral equation along L 


r 'dr ( x 1 ) 

7 l ■ _ g dx ' = ixl(x') (A3) 

where 


T indication that the principal part of the integral is to be 

J taken 

L union of smooth nonintersecting arcs or a group of dis- 

connected segments L' and L" along axis of reals 

l(x') prescribed boundary function along L that must be at least 
continuous within the open intervals L 1 and L" of L 
(l(x ) in the present case would correspond to prescription 
of v(x ) and v*(x’) along L' and L", respectively) 

) unknown function along L that is to be determined so as 

to be holomorphic in the region excluding L and approach 
on L the z(|) condition wherever prescribed (\|r(x') corre- 
sponds to v c (x') in the present instance) 

In reference 6, Muskhelishvili presents a thorough study of the 
integral equation (A5) , with particular attention given to the existence 
and uniqueness of solutions. The case where z(x') is real and imaginary 
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along alternating segments of L (mixed-type boundary values) has 
received detailed attention. It is this case that is of primary inter- 
est in the present analysis for the solution of v c along the given 

boundary and is now considered in detail. 

Specifically, the unique function v c (x',0 + ) = v(x',0 + ) + iv*(x',0 + ) 

satisfying equation ( A5) is needed; this function is to be holomorphic in 
the upper half plane y' >0, to be bounded at infinity, and to satisfy 
the given boundary conditions (see figs. 22 and 23 ) 

on L 1 , 


v(x',0 + ) = f(x') 


(A6a) 


on 


L" 


> 


v*(x',0 + ) = g(x') 


(A6b) 


where f(x') and g(x’) are continuous and finite along the open inter- 
vals of L’ and L" and are bounded and zero at one end point of each 
interval but may have -l/2 power singularities at the other end point of 
these intervals. In addition, it is assumed that g(x') for large values 
of x' is subject to the conditions 


lim g(x') 


x 


00 



g(x’) 


and 


g(x') - g(oo) < c ° - ^ . t . ai ^ (a > 0) (A7) 

(x') a 


where 


g(~) = lim g(x') 

X ’ ->+00 
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From the subsequent analysis it can be seen that inequality condition (A7) 
is satisfied in the present case. The general unique solution for v c 

satisfying conditions (A 6) and (A7) is given by Muskhelishvili (ref. 6 
p. 253) as * 


T.(x',ah =iE 

111 Jr . (x 1 ) 


U 


fiJD h(|) jR a (x') 

d| + Constant ■■ ■ — (a8) 




ftbU 7 ) 


where 


h(|) 


[f(5) on L' 
|ig(§) on L" 


= £ 0 - ) 


P 

R b (|) = U ( t - b.) 
j=i V jJ 


and a-jbj (j - 1; 2, . . . p) define the disconnected segments a b b b , 

a 2 b 2 , . . . along the axis of reals (see fig. 23). The application of 
the general solution (a 8) to the particular case herein requires some 
additional considerations. First note that in the present instance, 
along L ' , 


v(x’,0 + ) = f(x’,0 + ) = 0 (A9a) 

and along L", 


v*(x',0 + ) = g(x',Ot) =0 (A9b) 


Now consider a modification of the given boundary- value problem in which 
the singularity in v(x',0 + ) in the neighborhood above x ' is temporarily 
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neglected. For this case v approaches zero uniformly as x Q ' is 

approached from above. With this stipulation and conditions (A9) the 
general solution (A8) reduces to the complementary solution (or homo- 
geneous solution) 



(A10) 


where e^ is an arbitrary constant. The constant e^ has been evaluated 
in the body of the paper (see eq. (l6)). 

Consider now the contribution to the solution of the singularity 
in v (of the order -l) as Xq' is approached from above along x’ = x Q '. 

This contribution to the solution is represented by the integral term on 
the right-hand side of equation (A8) and can be evaluated by use of the 
following expediency. Assume the boundary is cut at x Q ' and then dis- 
placed an infinitesimal distance 2e parallel to itself; this open 
region is then filled in by a 2e segment of the x Q = x Q ' line that 

is normal to the boundary. Now the modified boundary in the infinitesimal 
neighborhood of x Q ' contains the singularity in v that was originally 

in the immediate neighborhood above Xq ' . Along this modified portion of 
the boundary between x Q ' - e and x Q ' + e the v-function takes on the 
form 


v(x 0 ' ± e) = 


Constant 
- a l)(“2 " 0 


(All) 


where 



and 


02 = 


+ e 


Equation (All) merely represents the transformations of the original -1 
singularity from the x' = x Q ' line to the y x = 0 line, the unit sin- 
gularity redistributing itself as two -l/2 power singularities, one at 
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x Q ' - e and one at x Q ' + e. Substituting the expression for v from 

equation (All) for h(| ) in equation (a 8) results in the following 
expression for the singularity contribution to the complex v-velocity: 


fv (x*)l = lim - 

1 

x' - a l x o lc k - b x 

Constant 


2 e— >0 

rti \ 

x ' -hJ ai=Xo .. £ 

_{(s - a i) (°2 - 0 

£ - x' 


(A12) 

Application of the mean- value theorem for the prescribed infinitesimal 
range of integration allows equation (A12) to be more simply expressed as 


llm . ~r 


0 rti , x . 


x' - a n 


na 2 =x 0 '+€ 


d£ 


a l =x o ' “ e ^ ~ " a l)( a 2 “ 0 


(A13) 


where ' is an arbitrary constant and the term 


'I - a. 


has been made 


1 


part of the constant. Integrating equation (A13) and then letting e >0 

yields 



(AlU) 


The constant e 2 has been evaluated in appendix B. 

This procedure for obtaining equation (Al4) for {v c (x'")] is by no 

means general; however, an evaluation of the integral of equation (A12) 
together with detailed considerations of the modification of the original 
boundary near x' = x Q ' (discussed previously) leads to the same expres- 
sion for F c (x'3 2 that is given by equation (Aik) . 

The complete solution for the complex v-velocity in the z' -plane is 
given by the sum of equations (A10) and (klk) . Transforming these rela- 
tions to the z -plane by lose of equation (A4) yields the following equation 
for the v-velocity in the z-plane along the boundary -x^ < x £ x^: 
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v c (x,0+ 



x ~ x h 
Xd - x 


(A15) 


The real and imaginary parts of equation (A15) correspond to equations (9) 
and (10) that are used in the determination of the expressions for the 
pressure velocity. 
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APPENDIX B 


EVALUATION OF CONSTANT e 2 


The constant e 2 may be evaluated by an integration of the w-velocity 
(the downwash velocity) along the rj-axis in the v-plane (fig. 5) from 
point d (the image of the wing leading edge) to unity (the end point of 
the Mach cone image) . Since the w-velocity is constant (equal to aV) 
along the wing segment 0 ^ t} ^ (d - e) £ and equal to zero at q = 1, 

then the integration of the downwash between the definite limits 

T) = (d - e) e — ^0 and- *1 = 1 is a known quantity and the constant & 2 _ 

is readily determined. 

The ensuing analysis makes use of the z-plane (fig. 6) which, as 
indicated in appendix A and the body of the paper, is obtained from the 
v-plane by a conformal transformation. In the z-plane the downwash 
velocity w(x,0 + ) in terms of the u-velocity may be expressed as follows 
(see eqs . (7) ) : 


u-velocity distribution over the wing-fin contour through the classic 
relation (from refs. 6 and 11) 


w(x,0 + ) = lim -B 






where 
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The functional notation — u(x,0 + ) is used to indicate 

v d y 


-d u(x,0 + ) = lim 
°y y— >0 4 


4 - u(x,y) 
oy 


Incorporating equation (B3) into 


equation (B2) results in the following expression: 


— u(x,0 + ) = - 

8y 


^ 1 + h 2 d 


( — 1 
l? 

x h + x ' 

Go 

jtB dx 


X - x' \ 

- x h 1 

x d - x ' 


dx ' + 


' x d 


- x h 


x' (x - x' ) 


1 


X h + X' 


dx’ 


x a 


- X' 


(-x h ^ X’ S XcLi 

Xd < x £ l) (B4) 


Evaluating the integrals in equation (B4) and differentiating the result 
with respect to x yields 


^1 + h 2 

x (3x h + Xd) - 2x h x d 


X + x h 

B 62 

2x2 ( x " x d) 

1 

X - x d 


— u(x,0 + ) = 
8y 


(x d < x < l) (B5) 

Insertion of equation (B5) for — u(x,0 + ) into equation (Bl) yields the 

°y 

following integral expression for the downwash: 


|w(x,0 + )J 


x \1 + h^ 

g e 2 

x d 2 


fnl f— 


ixd 


-(3 x h + x d) 


x h + x d 

XjX 


( x h + x d) ( x - x d 2 ) 2x h 


dx 


<iU) 


(B6) 


x d( x - x d) 
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where 


P(x) = (l - x 2 )^ - x d )(x - x h ) 

and f on the integral sign indicates that the finite part of the inte- 
gral must be taken. The evaluation of the integrals of equation (b 6) have 
been performed in appendix C. The following expression results for the 
downwash: 




+ E '( k ) 
k 


(B7) 


where K'(k) and E'(k) represent the complete elliptic integrals of 

the first and second kinds, respectively, each with modulus ]ll - k 2 . 
The quantity k is functionally defined as 


k = 


1 " x d x h - 


1 - x h 2 )( X ' x d 2 ) 


% " x h 


(B8) 


The downwash expression (B7) can be compactly expressed as 



(B9) 


Then eg is evaluated as 



- (-Va) = Va = 




1 + 


h 2 G 


e 2 


Va 1 


(BIO) 


40 


MCA TN 5O7I 


where G in the z-plane is defined by 


G = 


x h K'(k) + 



(Bll) 


In the original space, that is, the XYZ-space (see fig. l) , the expressions 
for k and G become 


k = 


1 - Bt 0 \|B 2 r Q 2 (l - B 2 t Q 2 ) + B^t" 2 - (l - B 2 t Q 2 )^ 
^B 2 r 0 2 (l - B 2 t Q 2 ) + B 2 t Q 2 - Bt 0 


1 - B 2 r D 2 


(B12) 


G = 


2k 


\ j B2 r 0 2 (l - B 2 t Q 2 ) 


+ B 2 t Q 2 - Bt 0 


Bt 0 K’(k) + 


E’(k) 


(B13) 


or 



(B14) 


(B15) 


The variation of k with BAy ^Ay is the aspect ratio of the vertical 
tail) for different values of BAg ( Ag is the aspect ratio of the 

horizontal tail) is presented in figure 7- Similarly, the variation 
of G with BAy for various values of BAjj may be obtained from fig- 
ures 8 and 9- 


NACA TN 3071 


Ul 


APPENDIX C 

EVALUATION OF INTEGRALS APPEARING IN THE EXPRESSION 
FOR THE DOWNWASH 




From appendix B the integral expression for the dovnwash 


is given 


[w(x,o + r) = — + h ~ 

x d 2 


(5 x h + x d) 


dx 


ix d \|(i - x2)(x - x d )(x - x h ) 


x h + x d 


dx 


x d Jx d X] j(l - X 2)(x - x d )(x - x h ) *** Jx d x 2 ^(l - x 2)( x - x d )(x . 


dx 


dx 


K + x d)(^ - x d 2 ) f p 

** Jx d (X - x d )3/2 p : x 2)( x . Xh) - 


(Cl) 


taLn^Th/* ^ ndiC f teS that the finite P ar t of the integral must be 
aken. The integrals appearing on the right-hand side of equation (Cl) 

ia + 7 el \ ipti< \ a ? d Can be transformed into the standard form 
(plus elementary integrals) by use of the linear transformation 


T = 


n k - ^ 

1 + r - X 

1 - toc^ 


k - 


Xd 


1 - toe. 


+ x 


The transforms for x d and x h are 


Xd = 


2k + k^x h + x h 
k 2 + 2kx h + 1 


Xv, = 


-2k + k 2 x^ + x^ 


^ - 2kx. 


d + 1 


(C2) 


(C3) 


J 
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where 


k 


1 


x h x a 




) 


x a " x h 


Substitution of equations (02) and (Cj) into the downwash expression, 
equation (Cl), yields 


[vt(t,0 + )J 


VI + h 


1/k 


2k j fa + * d )(l - k 2 ) 

- x h| k ( k - X d) 


dx 


l/k iJ(t 2 - l)(l - k 2 t 2 ) 


1 ( x h + x d)( X ~ ^ P T dt 

X d k5 J l/k - ^) 5/2 ^ 2 - 1 


K x h + x d)( L - k2 ) f 

nl 

dx 

x d k4 

1 — 1 

1 

CM 

CjcvT 

1 

CM 

X 

1— ) 


+ 


2 (i - x d 2 ) 
x d ( k2 * 2 x d k + x ) 



(C4) 


It is obvious from inspection that equation (C4) can be written in 
a more compact form; however, the presenteh form of equation (C4) is 
appropriate for rapia examination of limiting cases of the wing-fin 
arrangement (for example, when x h = 0, that is, when the fin aisappears, 

the corresponaing aownwash integral is reaaily ascertainea) . The inte- 
grals appearing in equation (C4) are reaaily evaluatea with the aia of 
the following substitutions: 


t = sn(u/k) 



cn(u/k) 




(C5) 


4 - r 2 k 2 = an(u/k) 
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where sn(u/k), cn(u/k), and dn(u/k) represent the Jacobian elliptic 
functions of argument u and modulus k. The evaluation of the inte- 
grals is as follows: 


TL 

P K 

dx 

du 

JiA \|t - ^ 2 )(* 2 - 1) J 

i 

K+iK' 


iK'(k) 

i 


-K'(k) 


pi 


u 


l/k 



dx 


"»K 

- dn 2 u du = E'(k) - K'(k) 

i 

JK+iK' 


In order to evaluate the finite parts of the integrals 



the following relation for the finite part was found most convenient: 


f pb 


U 


a 


A(x)dx 
(t 2 - a2) J / 2 


Pb 


A(x) - A(a) 

(t 2 - a 2 ) 3/2 


bA(a) 

a^b^ - 


(C 6) 


where A(x) is an integrable function in the closed interval ab. 
Studies of finite-part concepts are presented in references 15 and l6. 


44 


MCA TN 5071 


Application of equation (C6) results in the following evaluations of the 
preceding integrals: 


1 

b 3 

Tj 

1 — 1 

c 

«H 

f 1 dT k ^ | 

f 1 ) 

dr k 1 *- 

Jl/k ^2 . _ -l J 

iA (A - ^) 5 / 2 f ' 2 - 1 J 

1/k (x 2 k2 _ l) 3/2 ik ' 2 


= k 5 I JiL + iL - f K S^ du _ jL 
Jk+ik- u lk ' Jk+ik- ^ ik ' 2 


= ^§’00 - K'(kj| 


f ni 

T 

dr 

b 3 

r~ H 

r 

dx 

1 

~a 

dT 

k? 

U - $ 

1 

f— 1 

1 

C\J 

t- 

^ ( ,2 - $ 

i 3/2 \ 

I ’ 2 - 1 k, J 

l/k I 

^2 _ _L^ 
v k 2 y 

I OJ 
1 #r ^ 

(S_ 
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TABLE I.- SUMMARY OF PRESSURE COEFFICIENTS FOR VARIOUS LIMITING CASES 
OF THE GENERAL TAIL ARRANGEMENT SHOWN IN FIGURE 1 

Limiting 

conditions 

Pressure coefficient 

Remarks 


Slender tail 
(A h — ► 0; A V — » 0) 



11m G : 

Ay — * 0 

A H — >° 




/£p\ = M/i + 

W /v BG 




_^o I t Q + f2 + t 0 2 

+ *07 | \ r o 2 * t o 2 - Jr 2 + 1 


'&\ 0 _ 20 


(t), 


4*. ^ - t2 

\ <t „ 2 - tf | \jr 0 2 ♦ t „ 2 - ^r 2 ♦ 


(1 - ~ | ^ ~ ~ t2 

\ * t„ 2 ♦ 




Evaluated in 
ref. 12 


Vertical tall alone 
(A H = 0} 0 g BA V g 2) 


1 - fl - B2r Q J 


.fTS 


:’(k) 


:’(k) 


/^\ = 4p I r 

UJ V " BcVr - r 


Evaluated in 
ref. 13 


Leading edge of vertical tall sonic 
(M„ ,2; 0 f BA E S k) 




art .is 

Vl/V BO 


^(1 . B 2 ^ 2 ) . ^ 


Bt p + ^(l - B 2t 0 2) + B2t Q g 
1 - ^(l - B 2 ^ 2 ) + B2t 0 2 


V5£) „ . 2P 

BG 


Btp 




B2t 2 B2t2 


1 - B 2 ^ 2 1 - 


Btp 


B 2 tp 2 02 t 2 


^1 - B2tp2 1 1 - B2tp2 l - B2 t 2 


1 - B 2 ^ 2 |l - B 2 ! 2 1 - B 2 ! 2 


Btp 


1 Bt ° . 

J 

B 2 t2 

\,| f : bv 

h-^o 2 ■ 

L - B2 t 2 

ij, ■ ^tp 2 

J *V 

B 2 t 2 

y|l 1 - bV 

h - bV 

1 - B^t 2 


Leading edge of horizontal tail sonic 
(BA h .4; 0 g BA V < 2) 


lira C . 

BA H ->4 


or 

Wo r / 

/b 2 a v 2 | " 16 \ 

l 4 W ' J " 

4E-(Br 0 ) - 4E-^) 


-A 

U/V 0 


l 60 tp 


B^l - b\ 2 ^ 0 2 - r 2 


*a- 5 - 


1 - B2t2 


B\|l - B 2 t 0 2 | t2 ( 1 - B 2 ^ 2 ) + r p S 


Corresponds 
to preaaure 
loading on 
one panel of 
symmetrical 
triangular 
ving; evalu- 
ated in 
ref. 14 


Leading edges of horizontal and vertical tails sonic 
(BA V - 2; BA H = 4) 


11m G > 
BA V — >2 
BA«— +4 


(S* 75 *) 

f B 2 ky 2 J l6~ 


» Ka h 2 




I6p 




(n 




+ - B 2 ! 2 


Evaluated in 
refs. 14 
and 4 


Evaluated in 
ref. 4 


For the horizontal tail, pressure coefficients are presented only for the right panel. 

, l.fT^ 

V - ! *L_. 


BA V 
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TABLE II.- SUMMARY OF SIDESLIP DERIVATIVES FOR VARIOUS LIMITING 
CASES OF THE GENERAL TAIL ARRANGEMENT SHOWN IN FIGURE 1 


Derivative Remarks 


Slender tail 


(A H — 0; A V — > 0) 




Vertical tail alone 
(A h = 0; 0 < BA V ^ 2) 


nf2\ 

1 Ji bV 

Evaluated in 

1 - 1 1 - -r- 

BE'(k) 

ref. 13 

H5-J 



Cn fl = 
P 

5E'(k)BA v 





N)v - ' - 

K 4 

^BE'(k) 




— 


Leading edge of vertical tail sonic 
(BA y =2} 0 £ BA H g 4) 


0 B(4 + BA H ) 5 /2 


8)[2(3BA h + 4) 
110 3(4 + BA H )5/2 

(r \ -4\|2^ (trA" + 1 - HAt.) + ^ (\i 1 BAu^R" + sin“V^ ( 

( l * )v ■ 5ft ,fT^ (4 ; ^ [ (TOn n( n) H) 1* l- * ^ J; 

Leading edge of horizontal tail sonic 
(BA H =4; 0 £ BAy £ 2) 

«Ay 

Evaluated in 
ref. 14 


C = hn - 


/ \ 4Av 

(%)v = Tba^T 

V\ir) 




Leading edges of horizontal and vertical tails sonic 
(BA V = 2; BA h = 4) 

C I3 * - 1 

s*l (%)v=-fs W«-£i 

Can be evaluated 
from ref. 4 


-^5T 


H 



Figure 1.- Sketch of tail arrangement showing positive directions of 
velocities, forces, and moments. 
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Figure 2.- Orientation of tail arrangement for analysis. The horizontal 
surface, labeled wing, was originally the vertical tail. The vertical 
surface, labeled fin, was originally the horizontal tail. 
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Av(X,0±,Z) =0 
Au ( XjO^Z) =0 


u(X,Y,0 ± )=0 , 

v(X,Y,0 ± ) = 0 


v(X,0±Z) =0 
^»{X,0±,Z) = 0 
— u (X,0 ± ,Z) = 0 


w(X,Y,O i ) = aV 
|^u(X f Y,0±)-0 
^-v(X,Y,O i ) = 0 


u = 0 
v = 0 
w = 0 


Mach cone 


-Y,v 


u(X,Y,0 4 )*0 
v(X,Y,0±) = 0 


L-81263 

Figure 3 .- Sketch of wing-fin arrangement showing pertinent values of 
the Uj v } and w disturbance velocities and their spatial derivatives. 
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l-81261]. 

Figure 4.- Sketch of typical crossflow plane and its orientation in the 

X, Y, Z axes system. 
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a a ~o 
&/y *0 


Mach cone image 
_ u = ar=yu * o 


-/.0 



Mach cone image 

u^-nr~WjO 

10 



&U *0 

&OT--0 

1 -Lf 

Figure 5.- Sketch of wing-fin arrangement in the v = rj + i£ plane showing 
values of the u-, v- , and w- velocities and their derivatives along the 
real and imaginary axes . 
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L y 


VJl 

4=- 


-1.0 


Mach cone image dX. 

dyFx 


*-Left face 


Fm {upper half) 

-Right face — JD 


■+-Left face — >\*R/ghf face 
■< Fm joiner half) - 


h Q 


Mach cone image 
«===■ = = 

1.0 x 


(a) Transformed wing-fin arrangement in the z = x + iy plane. 
Numbers refer to numbers in figure 5* 


Mach cone image 
ns * O aj - o 


-Fm ( upper half) 
-Left face — *f <— Right face 


nr - o 


nr --0 




nr*-0 


Mach cone image 
0J--0 _ Jf =° 



Figure 6.- 


(b) Upper half of z-plane showing boundary values 
of v and v* along axis of reals. 

Sketch of z-plane including boundary values of v and v* 
along axis of reals. 
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Figure Variation of k with BA V for a range of BA H values. 
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values of BAg. 



B V 


16 

i B \ s 


values of 


- 1 with 
BA h . 


BAy 


for different 


Figure 9 >- Variation of 
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(a) Spanwise pressure distribution along a section k-m of vertical tail 



(b) Spanwise pressure distribution along a section l- k-1 of 

horizontal tail. 

Figure 10.- Spanwise pressure distributions along section k-m of vertical 
tail and section I-k-Z of horizontal tail. M = 1.17; A v =2; A H = 4. 


- 4.8 



8A V 

Figure 11.- Variation of BCy^ with BAy for various values of BAg. 
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Figure 12.- Variation of B^C^pjy with BAy 


for various values of BAjj. 
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Figure 13.- Variation of B (Cz p ) H with BA H for various values of BA y . 
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Figure lk.~ Variation of C n with BA y for various values of BA R . 
Moment about vertical axis through apex of tail. 
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Figure 15.- Variation of BCy^, B ( c ip) H > and with BA H 

for BA.y = 2 (leading edge of vertical tail is sonic). 



Figure 1 6 .- Variation of BCy^, B ( C *3)v' B ( Cz p) and Cn 3 Wlth Mv 
for BA h = ^ (leading edges of horizontal tail are sonic) . 
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Figure 17.- Illustrative variation of Cy^ with Mach number. A V = 2; 

Ajj = k. Results for supersonic leading edges obtained from reference k. 
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Figure 18.- 


A H 


= k. 


Illustrative variation of with Mach number. A V = 2; 

Results for supersonic leading edges obtained from reference 
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Figure 19.- Illustrative variation of C np with Mach number. A y = 2; 

A H = yawing-moment axis located a distance 2c r ahead of apex of 
tail. Results for supersonic leading edges obtained from reference 4. 
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Figure 20.- Plan-form modifications of basic tail surfaces for which the 
pressure expressions of the basic tail plan forms are also valid. 


/ 



Figure 21.- Wing — vertical-tail arrangements for which the pressure 
expressions of the basic tail plan forms are also valid. 
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Figure 22.- Sketch of the real axis of z '-plane shoving boundary values 


and singular points of v and v* along real axis, z' = 


2z - (x h + x d ) 



Figure 23.- Notation used in reference 6 in discussion of mixed- type 
boundary-value problem for the upper half plane. 
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